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ABSTRACT 


The  general  nonlinear  differential  equations  describing  the  interaction 
of  finitely  deformable,  polarizable,  heat  conducting  intrinsic  n-type  semi- 
conductors with  the  quasi-static  electric  field  are  transformed  from  the  un- 
known present  coordinate  description  to  the  known  reference  coordinate  descrip- 
tion, which  is  the  form  needed  in  the  treatment  of  problems.  From  the  dif- 
ferential form  of  each  balance  equation  in  the  reference  coordinate  descrip- 
tion, the  associated  integral  form  is  obtained.  The  resulting  integral  forms 
turn  out  as  expected  with  the  exception  of  the  one  due  to  the  balance  of 
linear  momentum  for  the  semiconducting  fluid,  in  which  an  important  change  in 
and  simplification  from  the  form  used  heretofore  is  introduced.  More  import- 
antly, the  previous  existing  integral  form  of  the  equation  of  the  balance  of 
energy  in  the  present  coordinate  description  is  transformed  to  a different 
form,  which  is  equivalent  to  the  original  form  only  when  the  field  variables 
are  differentiable.  The  revised  integral  form  in  the  present  coordinate 
description  is  then  transformed  to  the  reference  coordinate  description, 
from  which  an  energetic  jump  condition  across  a moving  non-material  surface 
of  discontinuity  is  obtained  which  is  consistent  with  all  the  other  jump 
conditions  obtained  from  the  other  integral  forms.  In  addition,  the  expres- 
sion for  the  quasi-static  electric  Poynting  vector  in  the  reference  coordinate 
description  is  determined. 


1.  Introduction 


In  a treatment  of  the  interaction  of  the  quasi-static  electric  field 
with  a finitely  deformable,  polarizable,  heat  conducting,  dielectric  continuum 
a system  of  rotationally  invariant  differential  equations  and  jump  conditions 
were  derived^.  First,  the  differential  equations  were  obtained  from  a 
systematic  application  of  the  laws  of  continuum  physics  to  a well-defined 
macroscopic  model.  Then,  the  integral  forms  obtained  from  the  respective 
differential  forms  were  taken  to  be  valid  across  moving  surfaces  of  discon- 
tinuity even  when  the  differential  forms  from  which  they  were  obtained  did 
not  hold  because  of  the  lack  of  differentiability  across  the  singular  surface. 
This  is  a standard,  perfectly  reasonable  procedure  for  obtaining  the  governing 
integral  forms  for  an  interacting  continuum.  The  jump  conditions  across  moving, 
not  necessarily  material,  surfaces  of  discontinuity  were  then  obtained  from 
the  respective  integral  forms.  The  resulting  integral  forms  and  associated 
jump  conditions  were  referred  to  the  unknown  present  (spatial)  coordinates 
of  material  points  whereas  it  would  have  been  more  useful  to  refer  them  to 
the  known  reference  (material)  coordinates  of  material  points.  However,  this 
is  not  a severe  limitation  since  it  is  a straight-forward  matter  to  transform 
from  present  to  reference  coordinates  of  material  points  using  well  known 
relations.  Nevertheless,  although  in  the  work^  in  question  the  appropriate 
reference  forms  of  the  jump  conditions  for  the  balance  of  mass,  linear 
momentum,  charge  and  entropy  may  readily  be  obtained  in  the  aforementioned 
manner,  the  reference  form  of  the  jump  condition  for  the  balance  of  energy 
may  not  be  obtained  directly  from  the  equations  presented  in  Ref.l.  Indeed, 
although  the  integral  form  of  the  equation  of  the  balance  of  energy  is  not 
incorrect,  the  present  form  of  the  jump  condition  for  the  balance  of  energy 
presented  there^  is  not  directly  useful  for  application  to  problems  concerning 


2. 


» 


the  propagation  of  surfaces  of  discontinuity  and,  hence,  is  not  actually 
correct.  This  dileitima  is  a consequence  of  the  fact  that  the  integral  form  of 
the  balance  of  energy  is  not  unique  because  the  differential  form  can  have 
terms  added  to  it  which  vanish  at  each  volumetric  point  and  yet  give  surface 
terms.  Since  the  field  varicJbles  are  not  differentiable  when  the  jump  condi- 
tions arc  obtained  from  the  integral  forms,  the  different  integral  forms 

2 

yield  different  jump  conditions  which  are  not  equivalent.  In  a later  work 

the  treatment  was  extended  from  the  case  of  deformable  insulators^  to  de- 

2 

formable  semiconductors.  In  the  latter  work  although  cin  integral  form  of 

the  equation  of  the  conservation  of  energy  is  written,  the  associated  jump 

condition  is  not  because  the  authors  were  aware  of  the  problem  involved. 

In  this  pap'er  the  proper  revised  differential  form  of  the  equation  of 

the  conservation  of  energy  is  found  for  the  more  general  case  of  the  deformable 

semiconductor.  From  the  revised  differential  form  the  appropriate  integral 

form  of  the  equation  of  the  conservation  of  energy  is  obtained  in  both  the 

present  £md  reference  coordinate  descriptions, from  which  the  required  correct 

energetic  jump  condition  in  the  reference  coordinate  description  is  obtained. 

One  by-product  of  this  work  is  the  determination  of  the  proper  expression  for 

the  quasi-static  Poynting  vector  in  the  reference  coordinate  description. 

In  addition  to  the  foregoing  a new  integral  form  of  the  equation  of  the 

conservation  of  linear  momentum  for  the  free-electronic  charge  continuum  is 

introduced.  The  new  form  is  better  for  obtaining  jump  conditions  because 

certain  assumptions  that  had  to  be  made  in  obtaining  jump  conditions  from  the 
2 3 

previous  form  ’ need  not  be  made  with  the  present  form.  Moreover,  the  3ump 
condition  at  a material  surface  of  discontinuity  resulting  from  the  present 
form  will  be  different  from  and  more  convincing  than  the  one  obtained  from 


the  previous  form. 


3. 


2 . Basic  Equations 

2 

It  has  been  shown  that  a system  of  general  nonlinear  rotationally  in- 
variant differential  equations  for  deformable  intrinsic  n-type  semiconductors 
can  be  written  in  the  form^ 


(T^.  + 

13 

ES  e 

T.,-p6..)  .=pdv, /dt  , 

13  ^ 13  ,1  ^3  ’ 

(2.1) 
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(v?  - ^ = p 0 dTl/dt  , 

(2.6) 

S ES 

where  t. . and  T. . denote  the  symmetric  parts  of  the  usual  mechanical  and 

13  13 

6 6 

free-space  Maxwell  electrostatic  stress  tensors,  respectively;  v^,  v,,  E^, 
^i’  ^i’  ^i  velocity  of  the  free  electronic  fluid,  velocity 

of  the  solid,  local  material  electric  field  exerted  on  the  free  electronic 
fluid,  Maxwell  electric  field,  electric  displacement,  electric  current  and 

6 6 6 IT 

heat  flux  vector,  respectively;  M'  , p , p,  cp,  cp  , ti,  [i  , 6 and  1]  denote  the 
free  electronic  charge  density,  free  electronic  pressure,  mass  density, 
electric  potential,  free  electronic  chemical  potential,  net  charge  density, 
residual  lattice  charge  density,  temperature  and  entropy  per  unit  mass, 
respectively;  and  is  the  permittivity  of  free  space.  The  motion  of  a 
point  of  the  solid  is  described  by  the  mapping 


yi=yi(X^t),  (2.7) 

where  y.  denotes  the  present  coordinates  of  material  points  and  X , the 
1 Xj 

reference  coordinates,  and  t denotes  the  time.  Clearly,  we  have 


4. 


V , = Sy . /9t  . 
1 1 


(2.8) 


In  (2.1)  - (2.8)  we  have  used  Cartesian  tensor  notation,  the  summation  con- 
vention for  repeated  tensor  indices,  the  convention  that  a comma  followed  by 
an  index  denotes  partial  differentiation  with  respect  to  a coordinate,  and 
the  convention  that  capital  indices  refer  to  reference  coordinates  and  lower 
case  indices  to  present  coordinates  of  material  points. 

The  associated  constitutive  equations  talce  the  form 


where 


and 


S SX  Sx 

ID  i,l^D,L  SE  ’ 1 o 1 aw 


'^i- ^1  = ^1,  A’  Ti=-ax/ae, 

P®  = (li.®)^ae®/9M.®  , 9®  = a (ti.®e®)/a)i®  , 


V = I <^1,  Al  ■ \ = ^i,  = ^i,L^i 


(2.9) 

(2.10) 

(2.11) 


(2.12) 


(2.13) 


In  addition,  we  have  the  conservation  of  mass,  which  may  be  written  in  the 
form 

PJ  = Pq  , (2.14) 

where  is  the  reference  mass  density  and 

J = det  yi  (2.15) 

At  this  point  we  note  that  the  above  system  of  equations  is  referred 
to  the  present  coordinates,  which  are  un)cnown,  eind,  as  a consequence,  it  is 
advantageous  to  refer  them  to  the  known  reference  coordinates.  To  this  end, 

Lj' 


analogous  to  the  Piola-Kirchhof f stress  tensor  F, 


which  is  defined  by 


5. 


’i  ■■■ 


n . T . . dS  = N . as  , 
1 13  L L3  o 


(2.16) 


4 5 

we  define  the  reference  electric  displacement  vector  3 by  ’ 

L 


n.D.  dS  = dS  , 

11  L L o ’ 


(2.17) 


where  dS^  and  denote  the  magnitude  of  and  unit  normal  to  an  element  of 


area  in  the  reference  configuration,  which  has  magnitude  dS  and  unit  normal  n. 


in  the  present  configuration.  By  virtue  of  the  well-]cnown  relation 


n . dS  = JX,  . N dS  , 
1 L,  1 L o ’ 


(2.18) 


from  (2.16)  and  (2.17),  in  the  usual  way,  we  find 


Lj  L,  1 13  ’ L L,  r 1 


(2.19) 


By  defining  the  reference  free-space  Maxwell  electrostatic  stress  tensor 


and  free  electronic  pressure  tensor  ^ . in  a similar  way  we  find 

hj 


ES  e 

M . = JX^  .T.  . , & . = .p 

t,3  L,ii3»  L3  L,3^ 


(2.20) 


Now, using  (2.19)  and  (2.20)  with  (2.1)  and  (2.2)2  and  employing  (2.14)  along 


7 

with  the  well-3aiown  identity 


(2.21) 


we  obtain 


(F  . + M.  . - ^ , = p dv . /dt , 
1*3  L3  L3  ,L  ^o  3 ’ 


= u.  . 
L,  L ^ ’ 


(2.22) 

(2.23) 


where 


p.  = Jp.  . 


(2.24) 


Equation  (2.22)  is  the  Piola-Kirchhof f form  of  the  stress  equations  of  motion 
for  the  conhined  deformable  semiconducting  continuum  and  Eq.  (2.23)  is  the 


I 


6. 


I 

f 

r 

I 

i 


reference  form  of  the  charge  equation  of  electrostatics.  From  (2.3),  the  j 

chain  rule  of  differentiation  and  (2.13)  , we  obtain  | 

•J  J 

w®  = (cp  + cp®)  (2.25) 

L,  , Li 


which  is  the  reference  form  of  the  equation  of  the  conservation  of  linear 
momentum  for  the  free  electronic  charge  continuum. 

In  order  to  refer  the  conservation  of  total  electric  charge  to  the  known 
reference  configuration  instead  of  the  unknown  present  configuration 
to  which  it  is  referred  in  (2.2)^^,  we  write 


the  substitution  of  which  in 
well-known  relations 


. / 


(2.2)^ 


ji  - , (2.  26) 

along  with  the  aid  of  (2.24)  and  the 


-1. 


d/dt  = (3/St)  + J '^k  k ~ dJ/dt 


(2.27) 


yields 


^ + J ^d|j,/dt  = 0 . 


(2. 28) 


Now,  we  define  the  reference  electric  current  vector  by 

V-  dS= 

which  with  (2.18)  yields 

k = ^^L,i^i- 

from  which  with  (2.28)  and  (2.21)  we  obtain 


(2.29) 


(2.30) 


(2.31) 


which  is  the  reference  form  of  the  conservation  of  total  electric  charge. 


7. 


Analogous  to  the  foregoing,  we  now  define  the  reference  heat  flux 


vector  Q by 

L 


which  with  (2.18)  yields 


n.q.  dS  = N Q dS  , 
1 1 L L O ^ 


= ^L.i'^i 


(2.32) 


(2.33) 


Substituting  from  (2.10)^^,  (2.13)  and  (2.33)  into  (2.6)  and  employing  (2.14) 

and  (2.21)  we  obtain 

A A O 

(2.34) 


-P^w®Q®  - = P^edn/dt, 


where 


p-  = Jp  . 


(2.35) 


Equation  (2.34)  is  the  reference  form  of  the  dissipation  equation.  The 
associated  reference  form  of  the  Clausius-Duhem  inequality  takes  the  form 

i * I 

where  F is  the  (positive)  rate  of  entropy  production. 

From  (2.4)^,  (2.9),  (2.10),  (2.19),  (2.26),  (2.30),  (2.33)  and 


(2.35)  the  pertinent  constitutive  equations  can  be  written  in  the  form 

^Lj  ^O^j,M  ;5E  ’ 0®L  " ^O  aw  ’ 

LM  L 


where 


5 = JX^  .E.  , 

L L,  1 1 ’ 


(2.37) 


(2.38) 


(2.29) 


and  the  constitutive  equations  in  (2.11)  and  (2.10)2  remain  unchanged  along 

6 G 

with  the  expressions  for  X,  e and  in  (2.12)  and,  of  course,  (2.5) 

and  (2.20)  are  to  be  employed.  From  (2.3)2  and  (2.13)2  we  note  that 


”k=-  '?,k- 


The  integral  forms  of  Eqs.  (2.22),  (2.23),  (2.25),  (2.31)  and  ( 

referred  to  the  known  reference  coordinates  of  material  points,  whi 
prime  interest  to  us  here,  clearly  may  be  written 

N (F  . + M . - ■^.)dS  = ^ I 0 V.  dV  , 

•J  L Ln  Xn  Lj  o dt  J o T o ’ 

S V 

o o 

J 

S V 

o o 

j w^dV^=J  N^(cptcp^)dS^, 


I N 2 dS  =-  — ndV  . 
J o dt  J ^ o ’ 

S V 

o o 


dV  i 0 . 
o 


The  integral  form  in  (2.43)  is  different  than  the  equivalent  equati 

2 3 8 

has  been  employed  heretofore  ’ ’ . Equation  (2.43)  is  a better  for 

obtaining  the  jump  condition  for  the  free-electronic  charge  continu 

both  material  and  nonmaterial  surfaces  of  discontinuity.  Since  Eq. 

0 

does  not  contain  the  free  electronic  charge  density  (i.  explicitly, 

0 

assumption  on  the  boundedness  of  (i.  need  be  made  in  obtaining  jump 
and,  as  a consequence,  the  resulting  jump  conditions  are  independen 
Maxwell  tensor  and  the  associated  assumptions  that  have  been  used  p 
From  the  integral  forms  in  (2.41)-  (2.45),  jump  conditions  across  m 
necessarily  material  surfaces  of  discontinuity  can  readily  be  obtai 


suitable  assumptions.  For  example,  at  a nonmaterial  surface  of  dis 


9. 


we  assume  that  all  volumetric  density  terms  remain  bounded  and  in  the  usual 
way  obtain 

N J =-  U„p^[vJ  , N r^_]  = 0,  [cp  + cp®]  = 0,  (2.46) 

\S^lI  = W®i  - Voi^I  ^ 

9 

where  is  the  intrinsic  velocity  of  the  singular  surface,  we  have  introduced 
the  conventional  notation  [a]  for  a - a^  and 

. + M . - -fff . . (2.48) 

L]  L]  L] 

Since  E.  remains  bounded,  we  have 

[cpl  =0,  (2.49) 

and,  hence,  from  (2.46)^  we  obtain 

jcp®]^  = 0.  (2.50) 

The  integral  form  of  the  equation  of  the  conservation  of  energy  and 
associated  jump  condition  has  not  been  discussed  here  because  the  existing 
integral  form  has  to  be  revised  in  accordance  with  the  detailed  discussion 
presented  in  the  next  section. 


3 . Integral  Fom  of  the  Balance  of  Energy 

In  this  section  we  consider  the  equation  of  the  conservation  of  energy 
for  the  deformable  intrinsic  n-type  semiconductor,  the  integral  form  of  which 
may  be  written 


d r ^ ^ e e 


1 ”1  ' i c e 

+ — s E,  E,  dv=  n.!T.,v,-p  V.  - 

2 o k kJ  J iL  1]  3 ^ 1 


h . - q.  - (v?  - V.  )|J.®e®  + V.  P,  E + V . i e E E,  i dS  . (3.1) 

11  1 1 ikk  i2okkJ 


where  denotes  the  electric  polarization  per  unit  volume  and  for  the  quasi- 


static electric  field 


10 


10. 


= cp[  Od  VSt)  + j^]  , 


(3.2) 


in  MKS  units.  Equation  (3.1)  is  the  extension  of  Eq.  (5.4)  of  Ref.l  for  the 
deformable  insulator  to  the  deformable  intrinsic  n-type  semiconductor  and  is 
also  Eq.  (4.14)  of  Ref . 2 when  reduced  to  the  special  case  of  the  intrinsic 
n-type  semiconductor  in  a quasi-static  electric  field.  Since  the  total  stress 
tensor  appears  in  the  equation  of  the  conservation  of  linear  momentum  for 
the  c<^mbined  continuum,  i.e.,  (2.1),  (2.22)  and  (2.41),  we  write 


where 


and  since 


S.  . = T . . + . - p®6.  . 

13  13  13  13 


h . = D.E.  - e E.E,  5..  , T.,=T^.-P.E. 

13  L 1 3 2 o k k 13J  ’ 13  13  13 


D.  = 6 E.  + P.  , 
1 01  1 


from  (3.3)-  (3.5)  and  (2.5)  we  note  that  in  addition 


^ S ES  e. 

S..  =T..  +T.,  -p6... 

13  13  13  13 


(3.3) 


(3.4) 


(3.5) 


(3.6) 


We  now  substitute  from  (3.3)  with  (3.4)^,  (2.4)^^,  (2.11),  (2.26),  (3.2)  and 

(3.5)  into  (3.1)  to  obtain 


d r ri 

. ' — Dv  V 

" - l2  ^ k 

V 

- D.E.v.  +E.D.V.  -cp“j'  - cp  -5-i  - cpj  . - q.  I dS 
133  331  1 ^ot  ^1  iJ 


6 G 1 

dt  . l2  + e + 2 Vk\. 

^ e SD. 


dV  = r n.[s. .V. 

■s  ^ 


(3.7) 


Now,  consider  the  term 


n.  [E  .D  . V.  - D.E  .V.  ] dS  = t , 

•11331  133  ’ 


(3.8) 


the  application  of  the  divergence  theorem  to  which  yields 

t = [E.D.v.  - D.E.V.]  , 

>'  331  133.1 


dV 


(3.9) 


11. 


Moreover,  appropriate  differentiation  and  the  substitution  of  (2.3)^  enables 


the  integrand  in  (3.9)  to  be  written 


(E.D.v.  -D.E.v.)  .=-  [cp(D.v.  -D.v.)  .]  ,, 


(D.v.  -D.v.)  ..  =0. 
31  1 3 ,13 


(3.10) 


(3.11) 


Substituting  from  (3.10)  into  (3.9),  employing  the  divergence  theorem  and 
(2.2)2,  we  obtain 

t = - n.cp(D.  .V.  + D.v.  . - D.v.  . - |a.v.)dS.  (3.12) 

J 1 1,3  3 1 3,3  D 1,  D ^ 1 


the  substitution  of  which  into  (3.7)  yields 


^ I pv  V + pe  + e E E 1 dv  = f n . fs  . . v . - cp(d . .v . 

it  J L2  k k 2 o k kJ  ^ iL  13  3 1, 3 ] 

SD.  V 

D.v.  .-D.v.  .-nv.  + -vr-  + 3.)  - 9 jf-q.  dS, 

1 3,3  3 1,3  1 1^  1 1-1 


(3.13) 


in  place  of  (3.7).  Equation  (3.13)  is  another  integral  form  of  the  equation 
of  the  balance  of  energy  referred  to  present  coordinates,  which  is  entirely 
equivalent  to  (3.7)  when  the  field  variables  are  differentiable.  However, 
when  the  field  variables  are  not  differentiable,  the  two  integral  forms 
obviously  are  not  equivalent. 

Equation  (3.13)  is  the  appropriate  form  for  obtaining  the  energetic 

jump  condition  in  the  reference  coordinate  description  while  Eq.  (3.7)  is  not. 

To  see  this  note  that  the  material  time  derivative  of  the  reference  electric 

displacement  vector  .5  defined  in  (2.19)  may  be  written 

L 2. 


dt 


where  (2.7)^  and  the  well-known  relation 


pBD.  -j 

. + v.D.  . - D.v.  . + V.  .D.  , 

, iL  at  3 1,  3 3 1,  3 3,3  1-1 


(3.14) 


(3.15) 


12. 


have  been  employed  in  obtaining  (3.14)  from  (2.19)2.  Now,  substituting  from 
(3.14),  (3.6),  (2.14),  (2.16),  (2.18),  (2.20),  (2.26),  (2.29),  (2.32),  (2.35) 

and  (2.48)  into  (3.13),  we  obtain 


I p V,  V,  + p e + e E,  E,  J I dv  = 1 N_  ( . v . 

dt  J L2  *^0  Ic  k ^o  2 o k k J o J lL  L]  3 


A * k)  - ’’*4  - »l]  <^0 . 


n: 


(3.16) 


which  is  the  integral  form  of  the  balance  of  energy  in  the  reference  coordi- 
nate description  we  have  been  after.  An  examination  of  (3.16)  reveals  that 
the  quasi-static  electric  Poynting  vector  referred  to  the  known  reference 
coordinates  takes  the  form 


h = ^ + 5 ) . 

L ^\dt  L ‘^L/  ’ 


(3.17) 


in  MKS  units.  The  energetic  jump  condition  across  moving  non-material 
surfaces  of  discontinuity  may  readily  be  obtained  from  (3.16)  by  assuming 
that  all  volumetric  density  terms  remain  bounded  with  the  result 


\yEj"j-  \ VoH- 


(3.18) 


where  h is  given  in  (3.17)  and 
L 


, e e , e 

.^1  u e 1 o 

Q = -r  V,  V,  + e + + - — E,  E, 

2kk  p 2pkk 


(3.19) 


It  should  be  noted  that  the  jump  condition  in  (3.18)  is  consistent  with  all 
the  other  j^p  conditions  appearing  in  Sec.  2,  while  the  jump  condition 
obtained  from  (3.1)  [or  (3.7)]  would  be  inconsistent  with  the  jump  condition 


in  (2.46)j^in  most  cases. 
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